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Starting from a quantum two-mode Bose-Hubbard Hamiltonian we determine the ground state 
properties, momentum distribution and dynamical evolution for a Bose Josephson junction realized 
by an ultracold Bose gas in a double- well trap. Varying the well asymmetry we identify Mott-like 
regions of parameters where number fluctuations are suppressed and the interference fringes in the 
momentum distribution are strongly reduced. We also show how Schroedinger cat states, realized 
from an initially phase coherent state by a sudden rise of the barrier among the two wells, will give 
rise to a destructive interference in the time-dependent momentum distribution. 
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Superconductor Josephson junctions are a paradig- 
matic example of macroscopic quantum coherence. The 
underlying physical mechanism is the Josephson ef- 
fect two superconductors connected by a weak link 
have coherent dynamical behavior determined by the rel- 
ative macroscopic phase of the superconducting conden- 
sates. Josephson junctions have been used to discuss 
fundamental concepts in quantum mechanics Q, per- 
form precision measurements and are now promising 
candidates to implement quantum information devices 
[3| • One important feature of superconducting Josephson 
junctions is the possibility to precisely control and adjust 
the state of the system by varying external parameters, 
eg a gate- voltage or magnetic flux. 

Bose Josephson junctions have been only recently pro- 
posed [H, and realized experimentally @, and many is- 
sues remain open. In the simplest configuration a Bose 
Josephson junction is realized by confining an ultracold 
Bose gas in a double-well potential. This configuration 
can be described using a two-mode model in which the 
bosons occupy the lowest level in each well. In the 
"classical" regime of large particle numbers and weak 
repulsive interactions the gas is well described by the 
Gross-Pitaevskii equation which, within the two-mode 
approximation, can be recast in the form of generalized 
Josephson equations for the time evolution of the relative 
phase and population imbalance among the two wells [j| . 
These equations differ from the original ones used for su- 
perconductor Josephson junctions [3| by the presence of 
a nonlinear coupling among the phase and population- 
imbalance variables. This term originates from boson- 
boson interactions in the mean-field approximation and 
gives rise to a rich dynamical behavior, displaying self 
trapping and 7r oscillations [5(. 

In this Letter we focus on the mesoscopic "quantum" 
regime beyond the Gross-Pitaevskii equation, in the limit 
of strong interactions and/or smaller values of N. This 
gets within reach of current experiments As inter- 
actions are increased phase fluctuations become more 



and more important while number fluctuations are sup- 
pressed; the ground state of the system approaches a 
regime which can be viewed as a mesoscopic Mott insula- 
tor. During the time evolution the phase coherence first 
degrades ("phase diffusion" [§]), but, in a closed quantum 
system, periodically revives, as it has been experimen- 
tally demonstrated [9]. At intermediate times between 
phase collapse and revival, Schroedinger cat states are 
predicted to form but are not easily observable in 
superconducting Josephson junctions 

The quantum behavior of superconducting Josephson 
junctions is usually accounted for by the standard phase 
model d, [IH . This model has been proposed to study 
the quantum fluctuations in a Bose Josephson junction 
[lij ]. and has been extended for large particle numbers 
N with subleading 1/N corrections lfj|; however it does 
not account for large population imbalance among the 
two sides of the junction. In this work we overcome this 
limitation. Using the quantum two-mode Bose-Hubbard 
Hamiltonian we investigate the ground state properties 
of the junction which we summarize in a "phase diagram" 
obtained by studying number fluctuations at varying well 
asymmetry and interaction strength. We characterize the 
phase-coherence of the junction by calculating the mo- 
mentum distribution. We then focus on the quantum 
dynamical behavior of the junction in a regime where 
Schroedinger cat states occur and show how they affect 
the time-dependent momentum distribution of the gas. 

Model We start by the two- mode Bose-Hubbard Hamil- 
tonian 
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where hi , a\ with i = 1 , 2 are bosonic field operators sat- 
isfying [hi, a}-] — Sij, Ef are the energies of the two wells, 
Ui > are the boson-boson repulsive interactions and K 
is the tunnel matrix clement, ic the Rabi oscillation en- 
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ergy in the case of a non-interacting model. The Heisen- 
berg equations of motion for this model yield 



ihd t ai = E®ai + Uifiihi — K&2 
ihd t ci2 = E%2 + U2n 2 ai — Kai, 
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where hi = ajfij. This is the quantum equivalent of the 
two-mode model used in the mean-field approximation 
ifli) = y/Niexp(i9i) and (hi) = Ni to describe Bosc- 
Josephson junctions 0. Indeed, by defining n = (Ni — 
N2)/2 and <fi = 62 — 6± the above equations are readily 
transformed into the Josephson-like equations 



hd t n = -2K\J (N/2) z - n 2 sin < 



hd t (j> = AE + nU s + Kn cos 0/ J (N/2) - n 2 (3) 



where AE = [E% + Ui(N- 1)/2- E$- U 2 (N- 1)/2] and 
U s = (U\ + 1/2). To go beyond the mean field model, we 
now transform the Hamiltonian (jTJ) into the exact quan- 
tum phase one by defining first the operators 17| hi 
\fn~i 
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e~ 1 \fh~i + 1 and then the relative- 
phase and relative- number operators e 1 ^ — £ ie2 e~ iei and 
h = (hi — h-x)j2. Throughout the paper we work at fixed 
total particle number N = hi + h% = N. In the new 
variables the Hamiltonian reads (up to a constant term) 
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with n = —AE/U S , related to the well asymmetry. 

Quasi- classical limit By using the commutation 
relations [y/N/2 + h + 1, e^] = (y/N/2 + h + l - 
y/ iV/2 + n)e^, and [\/N/2 - h + 1, eT^] 
(y/N/2 -h+l - y/N/2 - n)e-^, and by expand- 
ing the Hamiltonian |4]) for large TV we obtain: 

H = ^(n - n ) 2 -KN(l + ±- — ± 1/2 \ , os, 
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where cose 
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The leading term in this expansion corresponds to the 
standard phase model used for superconductors, which 
reads Hsjj = Ec(h — n g ) 2 /2 — Ejcos(f>, with Ec the 
charging energy, n g the dimensionless gate charge, and 
Ej the Josephson energy. Hence we have Ec — > U s and 
Ej — > ifiV. Also note that the well asymmetry no plays 
the role of n g , commonly used as an external param- 
eter to control the state of superconducting junctions. 
The subleading terms in Eq. (O correspond to a renor- 
malization of the plasma- Josepshon frequency \/U s KN, J y 



simplicity N even). By setting J x — (a\a,2 + a\ai)/2, 
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FIG. 1: Relative-number fluctuations (An } for the Bose- 
Josephson junction with N — 10 in the plane (720,7). 



such that in the limit U s — > tends to the Rabi fre- 
quency ojr. — Ej/N — K also present in the mean field 
solution. Eq. ||5j) shows how mesoscopic effects occur in 
Bose- Josephson junctions, involving a nonlinear coupling 
between h and <f>. 

Ground state of the quantum Hamiltonian. The quan- 
tum Hamiltonian ^ contains two very different regimes 
depending on the ratio 7 = KN/U S : for 7 3> 1 it yields 
a quasi-classical "superfluid" regime, where phase fluc- 
tuations are suppressed and the mean field approxima- 
tion applies, while for 7 <C 1 it yields the fully quantum 
" Mott-insulator like" regime where number fluctuations 
are suppressed (number-squeezed states). 

As we work in the two-mode approximation, the quan- 
tum Hamiltonian (j4| can be represented on a finite 
(N + 1) x (N + 1) matrix. In the Fock basis for relative- 
number states \j) the Hamiltonian (|3|) has a tridiago- 
nal form with matrix elements (j\H\j) = U s (j — no) 2 /2, 
(.] + l\H\j) = -Ky/N/2 + j +Ty/N /2-j and (j - 
l\H\j) = -Ky/N/2 -j + ly/N/2 + j. We have eval- 
uated the number fluctuations on the ground state of the 
system numerically, as represented in Figfl] in the plane 
(no 7 7)- The regions where number squeezing occurs are 
reminiscent of the Mott-insulator lobes of the phase di- 
agram of the Bose-Hubbard model 14|]; however we find 
only a smooth crossover between the two regimes since 
we are in the two-mode case. Note that for half-integer 
values of no the number squeezed regions are strongly 
suppressed even in the regime 7 <C 1. At these points 
the interaction energies of states j and j + 1 coincide 
favoring particle number fluctuations even if K is small. 

While in the limit 7 « 1 we have represented the 
Hamiltonian in the Fock basis |j) because its eigen- 
vectors are very close to the Fock states, in the oppo- 
site limit 7> 1, where large number fluctuations occur, 
it is more useful to represent the Hamiltonian by map- 
ping it onto angular momentum variables in the subspace 
at fixed J 2 = N/2(N/2 + 1) [H, EH (let us choose for 
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FIG. 2: Momentum distribution for a Bose Josephson 
junction for various values of the interaction strength 
(7 =100,1,0.1 and 0.001 from top to bottom), for TV = 10 
and no = 0. For the well wavefunction &o(x) we have chosen 
a gaussian profile of width a = O.ld. 



rewrite the Hamiltonian as: 

H = U s {J z -n f/2-2KJ x . 



(6) 



Its ground-state eigenvector is close to the angular- 
momentum coherent state UM 
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with a = tan(0/2) exp(— i<p) and J z \m) = m\m). Inter- 
estingly, the average energy on the state \a) is given by 



(a\H\a) = U s (l - l/N)n 2 /2 - 2Ky / (N/2) 2 - n 2 cos cj> 

(8) 

with n = —(N/2)cos8, which corresponds to 0(\/N) 
the mean field result for the energy, and at no = has a 
minimal value for 6 = 7r/2 and = 0. 

Momentum distribution The momentum distribution 
is one of the most accessible experimental observables. 
Here we consider the ground-state momentum distri- 
bution for the Bose- Josephson junction using the ex- 
act quantum phase model. The field operator in the 
two-mode approximation reads ^(x) = Y^i=i^i( x )^ii 
where $i(x) denotes the ground state wavefunction of 
the well i. The one-body density matrix, defined as 
pi(x,x') = (x)^ (x')) , where the average is intended 
over the quantum state of the system, for the two-mode 
model reads pi(x,x') = Y%j=i ®* 3 (x)$ l (x')(a]a l ). The 
momentum distribution is the Fourier transform with 
respect to the relative variable of the one-body den- 
sity matrix, n(p) = J dx J dx' exp(— «p(x — x'))pi(x,x') 
and becomes then n(p) — Y^ij=i ®j (p)^i(p)( a j a i) > wl th 
$i(p) being the Fourier transform of $i(x). For a sym- 
metric well or a weakly asymmetric situation we choose 



= $ (x - d/2), $ 2 (x) = <& Q (x + d/2), with d be- 
ing the interwell distance, and hence we obtain for the 
momentum distribution 



n(p) = \MpW(N + e-^(J+) + e*"<J_». (9) 

This is the generalization of the result derived by 
Pitaevskii and Stringari [la ] for a Bose-Josephson junc- 
tion using the standard phase model, and has also been 
used by Gati et al. 0] to quantify thermal decoher- 
ence in the experiment. In the quasi-classical regime 
7 > 1 we can evaluate the average in Eq.® using 
the coherent state ©, this yields n(p) = \&o(p)\ 2 {N + 
2y/ (N/2) 2 — n 2 cos(pd + 0); here we put n — such 
that n — for the ground state. Hence for 7> 1 we ex- 
pect interference fringes in momentum space fl5| while 
in the fully quantum regime 7 — > the matrix elements 
(J+,_) are vanishingly small, and the interferences are 
washed out as illustrated in FigJ^l Notice that this im- 
plies averaging over repeated measurements, as a single 
measurement would still yield interference fringes [2~q| . 

Schrddinger cat states Schrodinger cat states are quan- 
tum superpositions of macroscopic states. We suggest 
that such states might be realized as a result of the time 
evolution following a sudden rise of the barrier between 
the two wells, starting from an initially coherent state 
(ie in the regime 7 ^> 1). This procedure has the ad- 
vantage of starting from a quasi-classical state which is 
what is currently realized in experiments. For simplic- 
ity, we consider the symmetric case no = 0. If at time 
t = 0+ we set the inter-well coupling K in the Hamil- 
tonian to zero, then the time evolution is governed by 
the term U s J 2 /2 in the Hamiltonian. For each basis vec- 
tor \m) of the coherent state, the time evolution is given 
by \m(t)) = exp{-i2nm 2 t/T)\m), where T = Anh/U s is 
the revival period [2l[ such that \a(T)) — \a). Consider 
now the special times T/2q, q integer. In this case the 
phase factor governing the time evolution of the state \m) 
becomes exp(— mm 2 /q), which has the (anti-)periodicity 
property exp(—«7r(TO + g) 2 /(7) = (— l) 9 exp(— iirm 2 /q) de- 
pending on the parity of q. Hence for even q we can per- 
form a discrete Fourier transform to obtain a cat state 
given by a superposition of q coherent states, 
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where u k = (l/q)Yl 
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have |a(T/4)) = 
— a )]/y / 2, which is 



cat states exist for odd values of q. 

In particular, for q = 2 we 
[exp(-z7r/4)|a) + exp(i 7 r/4)(-l) Ar / 2 
a cat state given by the superposition of two coherent 
states. This kind of state was proposed for light coherent 
states by Yurke and Stoler [llj and for Josephson junc- 
tions by Gerry [l§]. However, in the latter case it was not 
easy to probe it. Here, we show that the cat states affect 
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for the unitary time evolution into the Schroedinger cat 
state is experimentally feasible as demonstrated already 
with the experiments in optical lattices 
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FIG. 3: Time evolution of (J x ) (a), and phase distribution 
at time t = T/8 (b) after a sudden quench of the coupling 
constants starting from an initially coherent state with a = 1 
for various values of N. 



the time-dependent momentum distribution. In partic- 
ular, by direct evaluation, we obtain that the contrast 
in the momentum distribution vanishes exactly for the 
two-component cat state. Furthermore, with increasing 
N, the contrast gets reduced on larger time regions as 
higher-order cat states develop at times different from 
t = T/4 (Fig|3Ja). For such cat states the phase distri- 
bution P(x) = ( a (0) exp(ix)\cx{t)) has equidistant equal 
peaks in the interval [0, 2tt] (Fig[3]b) which upon averag- 
ing strongly reduce the momentum-distribution contrast. 
We have also verified that these features are robust with 
respect to small tunneling among the two wells. 

Applications and experimental perspectives Number- 
squeezed states are particularly important for atom- 
optics applications, as their phase-diffusion time is longer 
than coherent states (ie usual Bose-condensates) Q . The 
measurement of the momentum distribution seems a 
promising route for the observation of quantum fluctu- 
ations and Schroedinger cat states on Bose Josephson 
junctions. One difficulty to overcome in the experimen- 
tal procedure is the precise control over the number of 
atoms in the junction, which should stay constant dur- 
ing the experiment. Atom losses may induce dephasing 
as discussed in 22]. On the other side, the time required 
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